----------------------------------------------------------------------------------------------------------------------------------
      name:  <unnamed>
       log:  C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_logs/appendix_table3.log
  log type:  text
 opened on:  24 Jun 2022, 10:08:50

.         
.         use "$Rep_smokelabor/1_build/regdata/county_quarter.dta", clear 

.         
.         local tbl_settings_log format(%6.3f) parentheses(stderr) asterisk()

.         local tbl_settings_pmil format(%6.1f) parentheses(stderr) asterisk()

.         local append replace

.         
.         reghdfe d_pc_qwi_payroll                        pm25 [aw=seer_pop] , a(fe_countyqtroy fe_styr) vce(cluster countyfip fe_
> stqtros)
(dropped 187 singleton observations)
(MWFE estimator converged in 8 iterations)

HDFE Linear regression                            Number of obs   =     74,725
Absorbing 2 HDFE groups                           F(   1,   1684) =      11.70
Statistics robust to heteroskedasticity           Prob > F        =     0.0006
                                                  R-squared       =     0.3215
                                                  Adj R-squared   =     0.2480
Number of clusters (countyfip) =      1,685       Within R-sq.    =     0.0016
Number of clusters (fe_stqtros) =      2,520      Root MSE        =   360.2342

               (Std. Err. adjusted for 1,685 clusters in countyfip fe_stqtros)
------------------------------------------------------------------------------
             |               Robust
d_pc_qwi_p~l |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
        pm25 |  -10.56581   3.089293    -3.42   0.001    -16.62507   -4.506551
       _cons |   249.7856   29.84766     8.37   0.000     191.2432     308.328
------------------------------------------------------------------------------

Absorbed degrees of freedom:
--------------------------------------------------------+
    Absorbed FE | Categories  - Redundant  = Num. Coefs |
----------------+---------------------------------------|
 fe_countyqtroy |      6676        6676           0    *|
        fe_styr |       631           0         631     |
--------------------------------------------------------+
* = FE nested within cluster; treated as redundant for DoF computation

.         summ pc_qwi_payroll [aw=seer_pop] if e(sample)

    Variable |     Obs      Weight        Mean   Std. Dev.       Min        Max
-------------+-----------------------------------------------------------------
pc_qwi_pay~l |  74,725  1.3907e+10    5687.648   3805.212   534.9928   62888.39

.         local ymean=`r(mean)'

.         regsave using "$Rep_smokelabor/2_analysis/output_tables/appendix_table3.dta", addlabel(yvar, "`ymean'") table(col_1, `tb
> l_settings_pmil') `append'
file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table3.dta saved

.         local append append

.         
.         reghdfe d_pc_qwi_payroll                        pm25 pm10 [aw=seer_pop] , a(fe_countyqtroy fe_styr) vce(cluster countyfi
> p fe_stqtros)
(dropped 312 singleton observations)
(MWFE estimator converged in 9 iterations)

HDFE Linear regression                            Number of obs   =     42,616
Absorbing 2 HDFE groups                           F(   2,   1070) =       8.40
Statistics robust to heteroskedasticity           Prob > F        =     0.0002
                                                  R-squared       =     0.3376
                                                  Adj R-squared   =     0.2537
Number of clusters (countyfip) =      1,071       Within R-sq.    =     0.0020
Number of clusters (fe_stqtros) =      2,519      Root MSE        =   376.3981

               (Std. Err. adjusted for 1,071 clusters in countyfip fe_stqtros)
------------------------------------------------------------------------------
             |               Robust
d_pc_qwi_p~l |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
        pm25 |  -13.64329   3.455862    -3.95   0.000    -20.42432   -6.862251
        pm10 |   1.853171   1.095907     1.69   0.091    -.2971992    4.003541
       _cons |   251.5041   38.13319     6.60   0.000     176.6798    326.3284
------------------------------------------------------------------------------

Absorbed degrees of freedom:
--------------------------------------------------------+
    Absorbed FE | Categories  - Redundant  = Num. Coefs |
----------------+---------------------------------------|
 fe_countyqtroy |      4162        4162           0    *|
        fe_styr |       631           0         631     |
--------------------------------------------------------+
* = FE nested within cluster; treated as redundant for DoF computation

.         summ pc_qwi_payroll [aw=seer_pop] if e(sample)

    Variable |     Obs      Weight        Mean   Std. Dev.       Min        Max
-------------+-----------------------------------------------------------------
pc_qwi_pay~l |  42,616  1.1635e+10    5975.228   3999.452   534.9928   62888.39

.         local ymean=`r(mean)'

.         regsave using "$Rep_smokelabor/2_analysis/output_tables/appendix_table3.dta", addlabel(yvar, "`ymean'") table(col_2, `tb
> l_settings_pmil') `append'
file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table3.dta saved

.         local append append

.         
.         reghdfe d_pc_qwi_payroll                        pm25 o3 [aw=seer_pop] , a(fe_countyqtroy fe_styr) vce(cluster countyfip 
> fe_stqtros)
(dropped 172 singleton observations)
(MWFE estimator converged in 9 iterations)

HDFE Linear regression                            Number of obs   =     64,248
Absorbing 2 HDFE groups                           F(   2,   1456) =       9.17
Statistics robust to heteroskedasticity           Prob > F        =     0.0001
                                                  R-squared       =     0.3248
                                                  Adj R-squared   =     0.2504
Number of clusters (countyfip) =      1,457       Within R-sq.    =     0.0018
Number of clusters (fe_stqtros) =      2,509      Root MSE        =   363.3500

               (Std. Err. adjusted for 1,457 clusters in countyfip fe_stqtros)
------------------------------------------------------------------------------
             |               Robust
d_pc_qwi_p~l |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
        pm25 |  -11.11598   3.377885    -3.29   0.001    -17.74202   -4.489936
          o3 |  -1.162037   1.063481    -1.09   0.275    -3.248155    .9240809
       _cons |   290.7749   33.91536     8.57   0.000     224.2468    357.3031
------------------------------------------------------------------------------

Absorbed degrees of freedom:
--------------------------------------------------------+
    Absorbed FE | Categories  - Redundant  = Num. Coefs |
----------------+---------------------------------------|
 fe_countyqtroy |      5748        5748           0    *|
        fe_styr |       631           0         631     |
--------------------------------------------------------+
* = FE nested within cluster; treated as redundant for DoF computation

.         summ pc_qwi_payroll [aw=seer_pop] if e(sample)

    Variable |     Obs      Weight        Mean   Std. Dev.       Min        Max
-------------+-----------------------------------------------------------------
pc_qwi_pay~l |  64,248  1.3383e+10    5763.877   3847.315   594.0308   62888.39

.         local ymean=`r(mean)'

.         regsave using "$Rep_smokelabor/2_analysis/output_tables/appendix_table3.dta", addlabel(yvar, "`ymean'") table(col_3, `tb
> l_settings_pmil') `append'
file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table3.dta saved

.         local append append

.         
.         reghdfe d_pc_qwi_payroll                        pm25 so2 [aw=seer_pop] , a(fe_countyqtroy fe_styr) vce(cluster countyfip
>  fe_stqtros)
(dropped 187 singleton observations)
(MWFE estimator converged in 10 iterations)
Warning: VCV matrix was non-positive semi-definite; adjustment from Cameron, Gelbach & Miller applied.

HDFE Linear regression                            Number of obs   =     40,363
Absorbing 2 HDFE groups                           F(   2,    978) =     863.30
Statistics robust to heteroskedasticity           Prob > F        =     0.0000
                                                  R-squared       =     0.3323
                                                  Adj R-squared   =     0.2485
Number of clusters (countyfip) =        979       Within R-sq.    =     0.0023
Number of clusters (fe_stqtros) =      2,502      Root MSE        =   380.3436

                 (Std. Err. adjusted for 979 clusters in countyfip fe_stqtros)
------------------------------------------------------------------------------
             |               Robust
d_pc_qwi_p~l |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
        pm25 |  -14.09433   4.381186    -3.22   0.001    -22.69193   -5.496719
         so2 |   -3.58851   2.329261    -1.54   0.124    -8.159435    .9824149
       _cons |   306.5599   39.98906     7.67   0.000     228.0857    385.0341
------------------------------------------------------------------------------

Absorbed degrees of freedom:
--------------------------------------------------------+
    Absorbed FE | Categories  - Redundant  = Num. Coefs |
----------------+---------------------------------------|
 fe_countyqtroy |      3873        3873           0    *|
        fe_styr |       628           0         628     |
--------------------------------------------------------+
* = FE nested within cluster; treated as redundant for DoF computation

.         summ pc_qwi_payroll [aw=seer_pop] if e(sample)

    Variable |     Obs      Weight        Mean   Std. Dev.       Min        Max
-------------+-----------------------------------------------------------------
pc_qwi_pay~l |  40,363  1.1322e+10    6114.536   4024.532   604.0206   62888.39

.         local ymean=`r(mean)'

.         regsave using "$Rep_smokelabor/2_analysis/output_tables/appendix_table3.dta", addlabel(yvar, "`ymean'") table(col_4, `tb
> l_settings_pmil') `append'
file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table3.dta saved

.         local append append

.         
.         reghdfe d_pc_qwi_payroll                        pm25 no2 [aw=seer_pop] , a(fe_countyqtroy fe_styr) vce(cluster countyfip
>  fe_stqtros)
(dropped 172 singleton observations)
(MWFE estimator converged in 10 iterations)

HDFE Linear regression                            Number of obs   =     31,534
Absorbing 2 HDFE groups                           F(   2,    784) =       4.51
Statistics robust to heteroskedasticity           Prob > F        =     0.0113
                                                  R-squared       =     0.3397
                                                  Adj R-squared   =     0.2524
Number of clusters (countyfip) =        785       Within R-sq.    =     0.0021
Number of clusters (fe_stqtros) =      2,487      Root MSE        =   388.9121

                 (Std. Err. adjusted for 785 clusters in countyfip fe_stqtros)
------------------------------------------------------------------------------
             |               Robust
d_pc_qwi_p~l |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
        pm25 |  -9.983638   4.006056    -2.49   0.013     -17.8475   -2.119771
         no2 |  -5.473712   3.945924    -1.39   0.166    -13.21954    2.272115
       _cons |   331.7128   57.81729     5.74   0.000     218.2178    445.2078
------------------------------------------------------------------------------

Absorbed degrees of freedom:
--------------------------------------------------------+
    Absorbed FE | Categories  - Redundant  = Num. Coefs |
----------------+---------------------------------------|
 fe_countyqtroy |      3054        3054           0    *|
        fe_styr |       625           0         625     |
--------------------------------------------------------+
* = FE nested within cluster; treated as redundant for DoF computation

.         summ pc_qwi_payroll [aw=seer_pop] if e(sample)

    Variable |     Obs      Weight        Mean   Std. Dev.       Min        Max
-------------+-----------------------------------------------------------------
pc_qwi_pay~l |  31,534  1.0549e+10    6211.774   4111.791   710.7908   62888.39

.         local ymean=`r(mean)'

.         regsave using "$Rep_smokelabor/2_analysis/output_tables/appendix_table3.dta", addlabel(yvar, "`ymean'")  table(col_5, `t
> bl_settings_pmil') `append'
file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table3.dta saved

.         local append append

.         
.         reghdfe d_pc_qwi_payroll                        pm25 pm10 o3 so2 no2  [aw=seer_pop] , a(fe_countyqtroy fe_styr) vce(clus
> ter countyfip fe_stqtros)
(dropped 127 singleton observations)
(MWFE estimator converged in 12 iterations)
Warning: VCV matrix was non-positive semi-definite; adjustment from Cameron, Gelbach & Miller applied.

HDFE Linear regression                            Number of obs   =     23,373
Absorbing 2 HDFE groups                           F(   5,    590) =     920.88
Statistics robust to heteroskedasticity           Prob > F        =     0.0000
                                                  R-squared       =     0.3460
                                                  Adj R-squared   =     0.2519
Number of clusters (countyfip) =        591       Within R-sq.    =     0.0034
Number of clusters (fe_stqtros) =      2,432      Root MSE        =   399.6524

                 (Std. Err. adjusted for 591 clusters in countyfip fe_stqtros)
------------------------------------------------------------------------------
             |               Robust
d_pc_qwi_p~l |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
        pm25 |  -15.39422   5.123587    -3.00   0.003    -25.45691   -5.331528
        pm10 |   2.592023   1.355567     1.91   0.056    -.0703005    5.254347
          o3 |  -2.460039   1.714901    -1.43   0.152    -5.828094    .9080151
         so2 |   1.210405   4.462758     0.27   0.786     -7.55442     9.97523
         no2 |  -7.587197   4.526828    -1.68   0.094    -16.47785    1.303461
       _cons |   431.1785   68.33794     6.31   0.000     296.9633    565.3938
------------------------------------------------------------------------------

Absorbed degrees of freedom:
--------------------------------------------------------+
    Absorbed FE | Categories  - Redundant  = Num. Coefs |
----------------+---------------------------------------|
 fe_countyqtroy |      2321        2321           0    *|
        fe_styr |       615           0         615     |
--------------------------------------------------------+
* = FE nested within cluster; treated as redundant for DoF computation

.         summ pc_qwi_payroll [aw=seer_pop] if e(sample)

    Variable |     Obs      Weight        Mean   Std. Dev.       Min        Max
-------------+-----------------------------------------------------------------
pc_qwi_pay~l |  23,373  9.4360e+09    6390.371   4244.521   710.7908   62888.39

.         local ymean=`r(mean)'

.         regsave using "$Rep_smokelabor/2_analysis/output_tables/appendix_table3.dta", addlabel(yvar, "`ymean'") table(col_6, `tb
> l_settings_pmil') `append'
file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table3.dta saved

.         local append append

.         
.         preserve 

.                 use "$Rep_smokelabor/2_analysis/output_tables/appendix_table3.dta", replace

.                 drop if var == "r2"
(1 observation deleted)

.                 drop if var == "_id"
(0 observations deleted)

.                 drop if strpos(var, "_cons_") > 0
(2 observations deleted)

.                 replace var="outcome_mean" if var=="yvar"
(1 real change made)

.                 replace var="observations" if var=="N"
(1 real change made)

.                 ingap 3

.                 ingap 6

.                 ingap 9

.                 ingap 12

.                 ingap 15

.                 list , sep(0)

     +-------------------------------------------------------------------------------+
     |          var      col_1      col_2      col_3      col_4     col_5      col_6 |
     |-------------------------------------------------------------------------------|
  1. |    pm25_coef   -10.6***   -13.6***   -11.1***   -14.1***   -10.0**   -15.4*** |
  2. |  pm25_stderr      (3.1)      (3.5)      (3.4)      (4.4)     (4.0)      (5.1) |
  3. |                                                                               |
  4. |    pm10_coef                  1.9*                                       2.6* |
  5. |  pm10_stderr                 (1.1)                                      (1.4) |
  6. |                                                                               |
  7. |      o3_coef                             -1.2                            -2.5 |
  8. |    o3_stderr                            (1.1)                           (1.7) |
  9. |                                                                               |
 10. |     so2_coef                                        -3.6                  1.2 |
 11. |   so2_stderr                                       (2.3)                (4.5) |
 12. |                                                                               |
 13. |     no2_coef                                                  -5.5      -7.6* |
 14. |   no2_stderr                                                 (3.9)      (4.5) |
 15. |                                                                               |
 16. | observations     74,725     42,616     64,248     40,363    31,534     23,373 |
 17. | outcome_mean     5687.6     5975.2     5763.9     6114.5    6211.8     6390.4 |
     +-------------------------------------------------------------------------------+

.                 saveold "$Rep_smokelabor/2_analysis/output_tables/appendix_table3.dta", replace
(saving in Stata 13 format)
(FYI, saveold has options version(12) and version(11) that write files in older Stata formats)
file C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_tables/appendix_table3.dta saved

.         restore 

. 
. log close
      name:  <unnamed>
       log:  C:\Users\ericzou\Dropbox\replicate_smokelabor/2_analysis/output_logs/appendix_table3.log
  log type:  text
 closed on:  24 Jun 2022, 10:08:54
----------------------------------------------------------------------------------------------------------------------------------
